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Broad cell tropism contributes to the pathogenesis of human cytomegalovirus (HCMV), but the 
extent to which cell type influences HCMV gene expression is unclear. A bespoke HCMV DNA 
microarray was used to monitor the transcriptome activity of the low passage Merlin strain of 
HCMV at 1 2, 24, 48 and 72 h post-infection, during a single round of replication in human fetal 
foreskin fibroblast cells (HFFF-2s), human retinal pigmented epithelial cells (RPE-1s) and human 
astrocytoma cells (U373MGs). In order to correlate transcriptome activity with concurrent 
biological responses, viral cytopathic effect, growth kinetics and genomic loads were examined in 
the three cell types. The temporal expression pattern of viral genes was broadly similar in HFFF-2s 
and RPE-1 s, but dramatically different in U373MGs. Of the 165 known HCMV protein-coding 
genes, 41 and 48 were differentially regulated in RPE-1 s and U373MGs, respectively, compared 
with HFFF-2s, and 22 of these were differentially regulated in both RPE-1 s and U373MGs. In 
RPE-1 s, all differentially regulated genes were downregulated, but, in U373MGs, some were 
down- and others upregulated. Differentially regulated genes were identified among the 
immediate-early, early, early late and true-late viral gene classes. Grouping of downregulated 
genes according to function at landmark stages of the replication cycle led to the identification of 
potential bottleneck stages (genome replication, virion assembly, and virion maturation and 
release) that may account for cell type-dependent viral growth kinetics. The possibility that cell 
type-specific differences in expressed cellular factors are responsible for modulation of viral gene 
expression is discussed. 
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Human cytomegalovirus (HCMV; species Human herpes- 
virus 5) is typically acquired in early childhood, although 
in developed countries this is often delayed. The virus 
establishes a life-long, persistent infection in myeloid pro- 
genitor cells (Kondo et al, 1994). The symptoms of pri- 
mary infection in healthy adults are generally mild, but 
infection during pregnancy can result in severe neurological 
damage to the fetus. HCMV is also a major contributor 
to morbidity and mortality among immunocompromised 



individuals. During the acute phase of primary infection, 
most body organs and a wide variety of cell types are infected 
(Sinzger et al., 1995). Indeed, the broad host cell tropism 
exhibited by HCMV in vivo is an important factor in its 
pathogenesis. However, the extent to which HCMV gene 
expression at the transcriptional level, and hence viral 
replication, is modulated because of the involvement of 
different cell types is unknown. In HCMV-permissive cell 
culture systems, both cytopathic effect (CPE) and infectious 
viral yield can vary markedly between cell types (Wang et al., 
2008). 
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The HCMV genome is organized into unique long and 
unique short regions (U L and U s ), each bounded by 
terminal and internal inverted repeats (TR l /IRl and TR S / 
IR S ). The genome of the low passage strain Merlin is 
235 646 bp in size and contains the estimated 165 protein- 
coding ORFs that are present in wild-type strains, although 
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two (RL13 and UL128) are interrupted by point mutations 
as a result of their instability in fibroblast cell culture 
(Dargan et at, 2010; Davison et al, 2003; Dolan et al, 2004; 
Stanton et al, 2010). Many HCMV genes are dispensable 
for growth in cell culture, and provide functions involved 
in processes such as cell tropism, immune evasion and viral 
temperance (suppression of viral growth) (Davison & 
Bhella, 2007). 

Regulation of HCMV gene expression is complex, and offers 
the potential for cell type-dependent regulation. Viral genes 
are transcribed by host RNA polymerase II in a regulated 
temporal cascade of immediate-early (IE), early (E), early 
late (E-L) and late (L) genes. Temporal regulation of 
some genes (e.g. UL4 and UL44) is mediated by alterna- 
tive promoters at different stages in the replication cycle 
(Isomura et al, 2007; Kerry et al, 1996, 1997). Transcripts 
that are antisense to viral mRNAs are produced in infected 
fibroblast and endothelial cells, and might have regulatory 
roles (Gatherer et al, 2011; Zhang et al, 2007). HCMV also 
encodes several miRNAs (Grey & Nelson, 2008), some of 
which may participate in control of viral or cellular gene 
expression (Stern-Ginossar et al, 2007). The HCMV gen- 
ome contains numerous binding sites for cellular transcrip- 
tion factors, and these may promote cell type-dependent 
differences in viral transcription kinetics and thus influence 
viral replication and pathogenesis. 

We investigated HCMV replication kinetics and synthesis of 
viral DNA and RNA in human fetal foreskin fibroblast cells 
(HFFF-2s), retinal pigmented epithelial cells (RPE-ls; Miceli 
et al, 1989) and astrocytoma cells (U373MGs; Kari et al, 
1992), as representative of cell types that are infected 
in vivo. Transcriptome profiling using DNA microarrays has 
historically provided a powerful technique for rapid analysis 
of global viral gene expression, and can help identify key 
regulatory effects above and beyond those due to gene- 
specific mechanisms. This technology has been employed for 
HCMV (Chambers et al, 1999; Yang et al, 2006) and several 
other herpesviruses (Aguilar et al, 2005; Ebrahimi et al, 
2003; Jenner et al, 2001; Kennedy et al, 2005; Li et al, 2006; 
Ohyashiki et al, 2005; Stingley et al, 2000; Tang et al, 2006; 
Tsao et al, 2009). Each of these studies investigated 
viral transcriptome activity in a single cell line, except for 
Kennedy et al. (2005), which found evidence for differential 
expression of varicella-zoster virus genes in two cell lines. 
We have developed a bespoke HCMV microarray platform 
to investigate transcriptome activity in three different cell 
types during a single round of replication by strain Merlin. 
We found that downregulation of certain virus genes may 
cause bottlenecks that operate at landmark stages in the 
HCMV replication cycle. 

RESULTS 

In order to correlate viral transcriptome activity in HFFF-2s, 
RPE-ls and U373MGs with concurrent biological responses, 
we examined viral CPE, growth kinetics and genomic loads. 



Viral growth kinetics 

The ability of strain Merlin to grow in HFFF-2s, RPE-ls 
and U373MGs after infection at high m.o.i. values was 
investigated. The extent of CPE at 72 h post-infection (p.i.) 
depended on the cell type (Fig. 1). All cells in the HFFF-2 
monolayers were rounded and clumped. In contrast, the 
U373MG and RPE-1 monolayers stayed intact, with cells 
remaining in contact with each other. Swollen and syncytial 
cells were a common feature of infected U373MG mono- 
layers, but were less frequent in RPE-1 monolayers, which 
exhibited little evidence of CPE. Thus, the cell type- 
dependent CPE exhibited by strain Merlin is similar to 
that reported previously for strains AD 169 and Towne in 
primary RPEs and astrocyte cultures (Detrick et al, 1996; 
Poland et al, 1990). 

Viral growth kinetics were compared in single-step growth 
experiments performed simultaneously with the three cell 
lines. The data were obtained from a single experiment, 
performed at the same time as the infections for microarray 
analysis (see below), using the same viral stock and batches 
of cells. Infectious viral yields from cells infected at an 
m.o.i. of 6 p.f.u. per cell were separated into cell-associated 
virus (CAV; Fig. 2a) and cell-released virus (CRV; Fig. 2b), 
and titrated on HFFF-2s. Compared with HFFF-2s, exit 
from the eclipse phase was delayed by 24 h in RPE-ls and 
U373MGs (Fig. 2a). Peak CAV titres were similar in HFFF- 
2s and RPE-ls, but were approximately 100-fold lower in 
U373MGs, and CAV accumulated most rapidly in RPE-ls 
and most slowly in U373MGs (Fig. 2a). CRV production 
was first detected at 72 h p.i. in HFFF-2s, but not until 
120 h p.i. in RPE-ls and U373MGs, and final yields were 
1000- and 10 000-fold lower, respectively, than in HFFF-2s 
(Fig. 2b). These results show that strain Merlin exhibits cell 
type-dependent growth properties similar to those of other 
strains (Detrick et al, 1996; Poland et al, 1990). 

To determine whether the three cell lines differ in 
susceptibility to viral infection, monolayers were infected 
with strain Merlin at 6 p.f.u. per cell and the numbers of 
infected cells were estimated by immunofiuorescent detec- 
tion of the UL44 protein. Compared with HFFF-2s (100 %), 
82% of RPE-ls and 77% of U373MGs were identified as 
infected. These small differences cannot by themselves 
account for the observed cell type- specific differences in 
viral growth. 

Viral DNA yield 

Whole cell DNA was harvested from HFFF-2s, RPE-ls and 
U373MGs either mock-infected (MI) or infected with 
HCMV at 6 p.f.u. per cell. Viral DNA yield was determined 
by real-time quantitative PCR (qPCR), and increased at a 
much slower rate in RPE-ls and U373MGs than HFFF-2s 
(Table 1). At 96 h p.i., when infectious viral yield ap- 
proached peak levels in each cell line (Fig. 2a), the viral 
DNA yield in HFFF-2s was 48- and 129-fold greater than 
in RPE-ls and U373MGs, respectively. It might be expected 
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Fig. 1. HCMV CPE at 72 h p.i. in HFFF-2s, 
RPE-1 s and U373MGs infected at 6 p.f.u. per 
cell. Ml, Mock-infected cells; HCMV, HCMV- 
infected cells. Bar, 500 \im. 



that the relatively low levels of HCMV DNA replication in 
RPE-ls and U373MGs (Table 1) contribute to delayed viral 
replication kinetics and lower viral yield (Fig. 2). However, a 
direct correlation cannot be drawn because the data in Table 
1 represent total viral DNA, whether packaged or unpack- 
aged. This suggested that cell type-dependent differences in 
viral maturation were likely to be important. 

Viral transcriptome activity 

A microarray platform specific for strain Merlin was 
developed to examine cell type- dependent differences in 
viral gene expression (i.e. transcript abundance) among 
HFFF-2s, RPE-ls and U373MGs. The data are presented in 
full in Table SI (available in JGV Online). Expression of 
approximately 96 % of the total number of viral ORFs was 
detected in HFFF-2s and RPE-ls, and 98% in U373MGs. 
Most ORFs were expressed at similar levels in HFFF-2s and 
RPE-ls (Fig. 3 and Table SI), with a temporally ordered 
increase in the numbers expressed. However, the pattern of 
expression in U373MGs was markedly different, with most 
ORFs, including those from genes exhibiting L kinetics in 



fibroblasts, expressed as early as 24 h p.i. (Table SI). Thus, 
compared with HFFF-2s and RPE-ls, the HCMV tran- 
scription cascade in U373MGs was either completed much 
earlier or significantly deregulated. 

Transcripts from a greater number of ORFs were detected at 
earlier time points in RPE-ls and U373MGs than in HFFF- 
2s, although those from a small number were first detected 
in HFFF-2s (Table 2 and Table SI). Totals of 41 and 48 
ORFs were differentially expressed (P<0.05) in RPE-ls and 
U373MGs (Table 2), respectively, compared with HFFF-2s. 
Of these, 13/41 ORFs in RPE-ls and 7/48 ORFs in U373MG 
were differentially expressed at multiple time points. In 
RPE-ls, all differentially expressed ORFs were expressed at 
lower levels than in HFFF-2s. However, in U373MGs, some 
ORFs were expressed at lower and others at higher levels 
compared with HFFF-2s, and a few exhibited a biphasic 
response (decreased at one time point and increased at 
another) (Table 2). Scrutiny of Table 2 reveals that 22 ORFs 
were differentially regulated in both RPE-ls and U373MGs. 

Differentially expressed genes were distributed among all 
temporal kinetic classes (Table 2). Twelve ORFs (UL44, 
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Fig. 2. Single-step growth curves of HCMV in HFFF-2s, RPE-1s 
and U373MGs infected at 6 p.f.u. per cell: (a) CAV and (b) CRV. 



Table 1. Estimates of HCMV genome loads 

Single experiment. Standard curve: slope=— 3.22; _R 2 =0.965. nd, Not 
done. 



Cell type 


Stage (h p.i.) 


Load (genomes per cell) 


HFFF-2 


MI 


<0.06 




24 


220 




48 


7 900 




72 


31000 




96 


62 000 


RPE-1 


MI 


<0.07 




24 


12 




48 


ND 




72 


270 




96 


1300 


U373MG 


MI 


<0.07 




24 


170 




48 


250 




72 


300 




96 


480 



UL80, UL85, UL86, UL89, UL94, UL97, UL98, UL99, 
UL104, UL114 and UL115) in RPE-ls and eight ORFs 
(UL50, UL55, UL69, UL85, UL98, UL99, UL114 and 
UL115) in U373MGs are responsible for conserved aspects 
of the replication cycle, in that they (core genes) were 
apparently inherited from the ancestor of all members of 
the family Herpesviridae. Among the families of HCMV 
genes (numbering 15) containing related members, the 
RL11 family (14 members) is strongly represented in the 
differentially expressed category: nine members in RPE-ls 
and five in U373MGs. In contrast, very few genes in the 
US22 family (12 members) were differentially expressed: 
one in RPE-ls and two in U373MGs. Five members of the 
US 12 family (10 members) exhibited enhanced expression 
in U373MGs. 

Concordance between microarray and Northern 
blotting data 

A selection of genes identified as differentially expressed in 
the microarray analysis were also assessed by Northern 
blotting (Fig. 4). 

IRS1 and TRS1 span the junctions between U s and TR S / 
IRs, and thus share identical promoters and 5 '-coding 
sequences, but have different 3 '-coding sequences (Chee 
et al, 1990). Two IRSl-related RNAs were capable of bind- 
ing the IRS1 microarray probe: a 3.5 kb full-length tran- 
script and a 1.7 kb transcript initiated within the ORF. The 
protein encoded by the latter has been shown to antagonize 
expression of that encoded by the former (Romanowski & 
Shenk, 1997). In Northern blots, levels of both transcripts 
were significantly reduced at all time points in RPE-ls 
compared with HFFF-2s and U373MGs (Fig. 4), thus 
confirming the results of the microarray analysis. 

UL4 produces a 1.8 kb transcript (Alderete et al, 2001). In 
Northern blots, this transcript was produced in reduced 
amounts in RPE-ls throughout infection and in U373MGs 
at 48 and 72 h p.i., compared with HFFF-2s, thus con- 
firming the results of the microarray analysis. 

UL82 and UL83 are expressed as a 3.5 kb bicistronic 
transcript terminating downstream from UL82 (Dal 
Monte et al, 1996). In Northern blots, the transcript level 
was similar in HFFF-2s and U373MGs to 48 h p.i., but 
moderately decreased (to 79 % of the level in HFFF-2s) at 
72 h p.i. in U373MGs. In RPE-ls, the transcript was less 
abundant throughout infection than in HFFF-2s, reaching 
54 % at 72 h p.i. These findings confirmed the results of the 
microarray analysis. 

UL93, UL94, UL95, UL96, UL97, UL98 and UL99 are 
expressed as a 3 '-co -terminal group of RNAs, which 
terminate downstream from UL99 and have L, L, E-L, E-L, 
E, E-L and L kinetics, respectively (Adam et al, 1995; 
Wing & Huang, 1995). In Northern blots, UL93 (10.5 kb 
transcript) and UL95 (7.3 kb) were expressed at similar 
levels in all three cell lines. UL94 (9.1 kb), UL96 (5.6 kb), 
UL97 (4.7 kb), UL98 (2.6 kb) and UL99 (1.6 and 1.3 kb) 
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Fig. 3. Scatter plots showing time-specific comparisons of HCMV 
gene expression in RPE-1 s and U373MGs relative to HFFF-2s. 
The solid line denotes equivalent expression levels and the broken 
lines indicate ± twofold differences, and dark and light grey points 
represent ORFs outside these boundaries. 



were expressed at lower levels at 48 h p.i. in RPE-ls than 
in HFFF-2s. UL96 (5.6 kb) and UL98 (2.6 kb) were ex- 
pressed at moderately higher levels in U373MGs than in 
HFFF-2s, and UL99 (1.6 and 1.3 kb) was produced in 
similar amounts in U373MG and HFFF-2s. These findings 
confirmed the results of the microarray analysis in regard 
to lower levels of expression of UL94, UL97, UL98 and 
UL99 in RPE-ls and higher expression of UL96 in 
U373MGs. However, results obtained for UL96 in RPE-ls 



and UL98 and UL99 in U373MGs differed from those 
obtained using the microarray. 

DISCUSSION 

Growth of HCMV strain Merlin was characterized in 
HFFF-2s, RPE-ls and U373MGs, which represent human 
fibroblast, epithelial and astrocyte cell types, respectively. 
Cell type-dependent differences were observed in viral 
CPE, infectious viral yields (CAV and CRV) and viral 
genomic yields. Viral growth in RPE-ls and U373MGs was 
characterized by delayed and reduced yields of CAV and 
CRV, relative to HFFF-2s (Fig. 2). Delayed viral growth 
kinetics have also been reported for primary RPEs infected 
with HCMV strain AD169 (Detrick et al, 1996). 

To investigate HCMV cell type-dependent transcriptional 
activity, a bespoke strain Merlin microarray was developed. 
Differential expression exhibited by 10 specimen genes 
was validated by Northern blotting. Minor discrepancies 
between the two approaches may be caused by several 
factors, including the higher sensitivity and statistical 
robustness of microarrays compared with Northern 
blots and the possible inability of microarray probes 
to discriminate fully among 3 '-co-terminal transcripts 
(Table 2; Chambers et al, 1999; Kennedy et al, 2005; 
Stingley et al, 2000). A total of 41 HCMV ORFs in RPE-ls 
and 48 ORFs in U373MGs were differentially expressed, 
compared with HFFF-2s. All differentially regulated ORFs 
in RPE-ls were expressed at lower levels, whereas some 
ORFs in U373MGs were expressed at higher levels and 
others at lower levels. The consequences of increased gene 
expression are unclear, but reduced expression of genes 
encoding regulatory or enzymic functions may play its part 
in restricting viral replication. 

Gross differences in viral entry, initiation of gene expres- 
sion or transition from E to L gene expression were ruled 
out as being responsible for differential expression of viral 
genes in RPE-ls and U373MGs, compared with HFFF-2s. 
The overall kinetics of gene expression were similar in 
RPE-ls and HFFF-2s, with expression of many ORFs being 
detected earlier in RPE-ls, and especially in U373MGs, 
compared with HFFF-2s (Table SI). Moreover, differential 
expression was not confined to genes of a particular 
temporal expression class. Despite viral DNA yields in 
RPE-ls and U373MGs being lower than in HFFF-2s, they 
appeared to be sufficient to maintain normal levels of 
transcription of most genes, including L genes. 

Functions associated with differentially 
expressed genes 

To shed light on the mechanisms by which viral growth 
in RPE-ls and U373MGs is suppressed, differentially 
expressed genes were grouped according to functional 
involvement at landmark stages in the viral replication 
cycle (Table 3). At each stage in RPE-ls and U373MGs, 
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Table 2. Differential expression of HCMV genes in RPE-1s and U373MGs compared with HFFF-2s 

Values are FCs in viral ORF signal intensity in RPE-ls/U373MGs relative to HFFF-2s: +, upregulated; — , downregulated. Only ORFs showing 
FCs>2 and P<0.05 are listed, nsd, No significant difference in signal between HFFF-2s and RPE-ls/U373MGs; SD, small difference in FC (<2). 
Selected 3'-co-terminal genes are in bold type. Protein properties and functions were adapted from Davison & Bhella (2007) by permission of 
Cambridge University Press. 



Gene 


Class 12 h p.i. 


24 h p.i. 


48 h p.i. 


72 h p.i. 


Protein properties and functions 


RPE-1 














RL5A 




* 


T 


/ .J 




ivi^ii idiiiiiy 


RL10 


E-L 


* 


* 


— 3.4 


NSD 


Virion pnvplnr>p (rlvrnnrotpin 


T?T 1 1 




* 


NSD 


Z 1 .7 


NSD 


RL11 family; IgG Fc-binding membrane glycoprotein 


RL12 








— 7 7 


NSD 


T? T 1 1 famin/' T*»ntcitn?p mpmhrati/ 1 cjiirmnrritpin 
IVij 1 1 idlllliy, L/LlLdllVC 111C111LI1 dllC LUU1U LC111 


RL13 


E-L 




NSD 


-17.3 


NSD 


RLll family, virion envelope glycoprotein 


UL4 


E-L 


-5.8 


-9.9 


-14 


NSD 


RL11 family; virion envelope glycoprotein 


UL5 


E-L 


NSD 


— 6.0 


—4.1 


NSD 


RT 1 1 Tamilv nnt^tivp fripirinranc nrotpi'ti 

Ivij 11 Idlllliy^ LJLlLdLlVC lllClllUldllC L/lt_/LClll 


UL10 


E-L 


* 


* 


-15.2 


NSD 


RLll family; putative membrane glycoprotein 


UL11 


E 




NSD 


—26.8 


NSD 


T?T 1 1 Tamilv mpmnrjinp cnvrnnrntpin 

IvL 11 Idlllliy^ 111C111L/1 dllC iLl VLUU1UIC111 


UL17 


E 


NSD 


NSD 


— 24.6 


NSD 




UT 22A 




t 


NSD 


— 8.2 


NSD 


Cpf i-ptp/l T? ANTPS-nitiHiricr nmtpin 

v.'LLl L KU 1VZ1.1N i IjU Ulim.llliL L/ltJLClll 


TTT 9^ 

U.LZ J 


L 




■f 

T 


NSD 


J. y 


UL25 family; tegument phosphoprotein 


UL41A 




t 




— 3 ? 
J.Z 


NSD 


i ULdLlvc lllcIllUI dllc Dl (JLclll 


T TT AA 


E-L 


t 


NSD 


Z.O 


NSD 


Processivity subunit of DNA polymerase 


T TT AQ 


E-L 


NSD 


/ .J 


NSD 


L.y 




TTT RD/TTT SO ^ 


L 




4- 


— 14 1 




r IUlcd.sc dllU CdDMU aCdlKJlU. UIULclIlt* 


UL82 


E-L 


NSD 


— 2.9 


— 15.7 


— 5.7 


1~)T TT? P tciTYiilv t Aon m An t nriAcnnfiTif Atpin nn71 ■ 

LJ VJ JTVJT ld.llll±y3 LCtUlllCllL L71HJ3 IJ11U L>1 U LC111 UU; 1 } 

transcriptional activator; targeted to ND10 domains; 
targets Rb proteins for proteosomal degradation 


UL83 


E-L 


t 


t 
\ 


—41.6 


— 5.4 


1~)T TT? P finnilv tponmpnt nrinQnriftnrntpin nnftS - QiinnrpQ'jp^ 

i~s LJ Ivl Ldlll.il LCilLLlllClll L'lHJoLJllV_JL'lV_JLClll LJL/U—/, □LiUUlCaoto 

111 Lei lei Ull I CapUllac 


UL85 


E-L 




± 

T 


— 19.6 


NSD 


("'ninririnpnt nf intprratriAmprir trinl pyps in ranQiH'; 

\_j\_JlllL/t_/llCll L VJL UlLClCdL'jUlllCllC Lll L/ltACa 111 CdLJcjlVJ-j 


UL86 


E-L 




* 


— 14.1 




A/Tuir^t" rancin nrntpiri ■ rrdnnnTiPTit nf tiPYnnc ann npntr^nc 

IVldlUl L.d.L'SlU. IJlULClllj U(J111IJ(J11C11L \JL 11CA.U11S allU IJCllLUllo 




E-L 


* 




— 5 0 




i LlLd-Llvc nil dsc aU.UU.lllL (Jl Lei lllllldac 


UL91 


L 


t 


NSD 


-4.9 


NSD 




UL94 


L 


* 




-14.7 


NSD 


Tegument protein; binds ssDNA 


UL97 


E-L 


t 


NSD 


—4.2 


NSD 


Serine-threonine protein kinase; phosphorylates UL44 


UL98 


E-L 


NSD 


— ins 


— 12 1 


_ 3 Q 
J. 7 


DlN'ase 


TTT QQ 


L 


NSD 




— Q Q 
O.O 


J.J 


IVlyi loLyldLcU. LcgUlllCllL Ul (JLclll pUZO 


UL104 


E 




* 


± 




\_-dpt>lU. UUl Ldl UlULclll 


UL1 12 


E 


NSD 


NSD 


—4.3 


NSD 


(~lt*f"npCT"i"ai"AC TlT\T A rptiliratiAri TM*r*T"piTic at \nt"dl i*Anli/~atir\ri 

\_/1C11Cc>L1 dlCS lylN Ji 1 C IJllCdLlUll L7HJLC111S dl Vlld.1 1 C IJllCdLlOll 

compartments 


UL114 


E 


*■ 


T 




NSD 


T Tt*o r*i 1 Tl NT A rrl"\^/~/"* citt ocp 

uidcii i_ji>x\ yiycutiyidac 


UL115 


L 


t 


-2.0 


-3.9 


NSD 


Virion envelope glycoprotein gL 


UL117 


E 


t 


NSD 


— 9.6 


NSD 


T?pniiirpn fnr virsl rpnliratinn rnmnartmpnts 

IVtU LL11CU- 1U1 Vll dl 1C L'llCd L1U11 CVJ111 Ual LlllCllLo 


UL1 19 


E-L 




NSD 


— 3.0 


NSD 


T<rf^ T<f"- ni n/ni net m atyi nran a (rlAr/~r\nt"r\tAin rplatpn trt C~\^C-') 

liiVJ J. C UlllUlllli 111C111U1 dllC iiiy COIJl OLClll ICldlCU LU VJA i. 


UL122 


TT-? T 


4. 

T 


NSD 


-36.6 


-6.1 


IE transcriptional activator IE2; interacts with basal 
transcriptional machinery and cellular transcription 
factors; specific DNA-binding protein 


UL132 




NSD 


-3.3 


-4.4 


NSD 


Virion glycoprotein 


UL148 




-3.6 


NSD 


-4.5 


NSD 


Putative membrane glycoprotein 


UL147A 




-3.1 


NSD 


-10.8 


-4.2 


Putative membrane protein 


UL146 






* 


* 


-34.4 


CXCL family; putative secreted CXC chemokine 


UL144 




t 


NSD 


-8.8 


NSD 


Putative membrane glycoprotein; TNF-receptor 
homologue 


UL141 






NSD 


-3.8 


NSD 


UL14 family; membrane glycoprotein; inhibits NK cell 
cytotoxicity by downregulating CD 155 


UL139 




* 




t 


-57.1 


Putative membrane glycoprotein 


IRS1 


IE 


-4.2 


-7.4 


NSD 


NSD 


US22 family, IE transcriptional activator, tegument 
protein, involved in host cell shut-off 
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Table 2. cont. 



Gene 


Class 


12 h p.i. 


24 h p.i. 


48 h p.i. 


72 h p.i. 


Protein properties and functions 


U373MG 














RL11 


E-L 


t 


4.9 


NSD 


SD 


RL11 family; IgG Fc-binding membrane glycoprotein 


RL13 


E-L 


t 


2.2 


NSD 


NSD 


nT 1 1 f "1 1 1 . • 

RL11 family, virion envelope glycoprotein 


UL4 


E 


NSD 


-3.1 


-3.6 


NSD 


RL11 family; virion envelope glycoprotein 


UL10 




t 


t 


2.5 


NSD 


RL11 family; putative membrane glycoprotein 


UL11 


E 


NSD 


NSD 


3.9 


NSD 


RL11 family; membrane glycoprotein 


UL13 


E 


NSD 


NSD 


2.8 


NSD 


Putative secreted protein 


UL19 




t 


t 


18.7 


NSD 




UL21A 




t 


3.1 


44.1 


NSD 


Facilitates viral DNA synthesis and late accumulation of 
IE transcripts 


UL22A 




t 


NSD 


—2.8 


NSD 


Secreted RANlES-binding protein 


UL25 


E 


t 


t 


4 


-4.6 


UL25 family; tegument phosphoprotein 


UL26 




t 


t 


21.9 


NSD 


US22 family; tegument transactivator of MIEP 


UL32 


E 


t 


t 


NSD 


-4.1 


Major tegument phosphoglycoprotein ppl50; binds to capsids 


UL33 


E 


t 


t 


25.1 


NSD 


GCPR family; virion protein; putative chemokine receptor 


UL34 




t 


NSD 


5.7 


NSD 


Repressor of US3 transcription 


UL38 


IE-E 


t 


t 


6.1 


NSD 


Anti-apoptotic protein; protects cells from ER stress 
induced by unfolded protein response 


UL43 


L 


t 


NSD 


3.7 


NSD 


US22 family; tegument protein 


UL49 


E-L 


NSD 


NSD 


NSD 


—2.2 




UL50 




t 


2.5 


7.5 


NSD 


Inner nuclear membrane protein; facilitates capsid egress 
from nucleus 


UL55 


E-L 


t 


3.3 


NSD 


NSD 


Viral envelope glycoprotein gB; mediates membrane 
fusion during entry 


UL69 


E-L 


t 


NSD 


5.4 


NSD 


Regulatory protein; tegument protein; contributes to cell 
cycle block; involved in mRNA export; exhibits 
nucleocytoplasmic shuttling 


UL82 


E-L 


NSD 


NSD 


NSD 


-8.8 


DURP family; tegument phosphoprotein pp71; 
transcriptional activator; targeted to ND10 domains; 
targets Rb proteins for proteosomal degradation 


UL83 


E-L 


t 


t 


NSD 


— 8.4 


DURP family; tegument phosphoprotein pp65; suppresses 
interferon response 


UL85 


E-L 


t 


NSD 


-2.2 


NSD 


Component of intercapsomeric triplexes in capsids 


UL96 


E-L 


t 


NSD 


2.8 


NSD 


Tegument protein 


UL98 


E-L 


NSD 


NSD 


NSD 


4.7 


DNase 


UL99 


L 


NSD 


NSD 


NSD 


3.8 


Myristylated tegument protein pp28 


UL112 


E 


NSD 


4.2 


3 


NSD 


Orchestrates DNA replication proteins at viral replication 
compartments 


UL114 


E 


t 


t 


2.8 


NSD 


Uracil-DNA glycosylase 


UL115 


L 


t 


2.1 


NSD 


NSD 


Virion envelope glycoprotein gL 


UL116 


E-L 


t 


NSD 


3.3 


NSD 


Putative membrane glycoprotein 


UL119 


E-L 


t 


NSD 


4 


NSD 


IgG Fc-binding membrane glycoprotein related to OX-2 


UL122 


IE-L 


t 


NSD 


NSD 


-5.6 


IE transcriptional activator IE2; interacts with basal 

' . ' 1 1* 1111 • • 

transcriptional machinery and cellular transcription 
factors; specific DNA-binding protein 


UL128 


E-L 


t 


t 


13.4 


NSD 


Viral envelope protein; cell tropism factor for endothelial, 
epithelial and leukocyte cells 


UL148 




NSD 


2.8 


NSD 


NSD 


Putative membrane glycoprotein 


UL147A 




NSD 


5.3 


3.4 


-3.8 


Putative membrane protein 


UL147 




NSD 


6.7 


8.1 


-2.6 


CXCL family; putative secreted CXC chemokine 


UL146 




t 


t 


t 


-63.7 


CXCL family; putative secreted CXC chemokine 


UL145 




t 


t 


13.1 


NSD 


RL1 family 


UL140 




t 


t 


NSD 


-5.3 


Putative membrane protein 


UL139 




t 


t 


NSD 


-15.5 


Putative membrane glycoprotein 


UL138 




NSD 


NSD 


2.6 


NSD 


Regulator of viral latency 


UL136 




t 


NSD 


3.5 


NSD 


Putative membrane protein 
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Gene 


Class 


12 h p.i. 


24 h p.i. 


48 h p.i. 


72 h p.i. 


Protein properties and functions 


UL135 




NSD 


NSD 


4.3 


NSD 


Putative secreted protein 


US12 


E 


NSD 


NSD 


4.3 


NSD 


US12 family; putative multiple transmembrane protein 


US13 


E 


NSD 


NSD 


6.4 


NSD 


US12 family; putative multiple transmembrane protein 


US16 


E 


t 


NSD 


8.6 


NSD 


US12 family; putative multiple transmembrane protein 


US19 


E 


NSD 


NSD 


5.7 


NSD 


US 12 family; putative multiple transmembrane protein 


US20 


E 


NSD 


2.9 


3.9 


NSD 


US 12 family; putative multiple transmembrane protein 



*Signal below detectable limits in HFFF-2s, RPE-1 and U373MGs. 
fSignal not detected in HFFF-2s but detected in RPE-ls/U373MGs. 
tSignal detected in HFFF-2s but not RPE-ls/U373MGs. 



expression of multiple viral genes was affected, thus 
demonstrating the broad effects of the underlying processes 
(Table 2). Downregulation of certain herpesvirus core genes 
in RPE-ls is compatible with a trend to temperance of virus 
replication. The situation is less clear for U373MGs, as 
differences in the levels of gene expression were generally 
smaller and many ORFs were upregulated, making the 
consequences for viral replication less predictable. However, 
most genes downregulated in U373MGs have roles in viral 
maturation and release. Many more replication steps were 
potentially impaired in RPE-ls than in U373MGs (Table 3), 
and yet higher infectious viral yields were generated (Fig. 2), 
implying a secondary obstacle to viral replication in 
U373MGs that has a greater impact on viral growth than 
impaired maturation. The apparent deregulation of the viral 
transcription cascade may be an important factor, likely to 
perturb post-transcripional events, including translation, 



post-translational processing, intracellular translocation of 
viral proteins and viral DNA synthesis, and thus contrib- 
uting to the low infectious viral yields. In summary, a series 
of potential bottleneck stages (Table 3) that may impede, but 
not halt, progression through the viral replication cycle was 
proposed in RPE-ls and U373MGs. This may account for 
the characteristic cell type-dependent viral growth kinetics 
observed. 

HCMV transcriptome expression was investigated over the 
period 12-72 h p.i., during which there was a 24 h delay in 
production of CAV in both RPE-ls and U373MGs, 
compared with HFFF-2s. In RPE-ls, and to a lesser extent 
in U373MGs, CAV yield increased at 72 h and 96 h p.i. (Fig. 
2a). It is possible that, as the protein products of under- 
expressed genes reach a critical level, replication bottlenecks 
are cleared sequentially. Moreover, accumulation of CAV in 



HFFF-2 RPE-1 U373MG 

12 24 48 72 12 24 48 72 12 24 48 72 



(a) 



-3.5 
•1.7 



(b) 
(c) 



1.8 



-3.5 




(e) 



f« M GAPDH 



Fig. 4. Northern blots showing representative 
HCMV transcripts produced in HFFF-2s, RPE- 
ls and U373MGs at 12, 24, 48 and 72 h p.i. 
(a) IRS1 (3.5 and 1 .7 kb), (b) UL4 (1 .8 kb), (c) 
UL83 (3.5 kb), (d) UL93 (10.5 kb), UL94 
(9.1 kb), UL95 (7.3 kb), UL96 (5.6 kb), UL97 
(4.7 kb), UL98 (2.6 kb) and UL99 (1 .6 and 
1 .3 kb) and (e) glyceraldehyde 3-phosphate 
dehydrogenase (GAPDH). Genome co-ordin- 
ates (GenBank accession no. AY446894.2) 
of the HCMV probe sequences were: IRS1 , 
198 271-197 985; UL4, 13 936-142 429; 
UL83, 1 21 565-1 21 800; and UL99 1 45 796- 
146 085. 
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Table 3. Differentially expressed HCMV genes grouped according to functional class in the replication cycle 



Expression is given in relation to HFFF-2s. No gene in these classes was overexpressed in RPE-ls. 



Functional class RPE-1 underexpressed U373MG overexpressed U373MG underexpressed 


Regulation of gene expression IRS1 


UL26 


UL122 


UL122 


UL69 


UL82 


UL82 






Development of replication compartments UL112 


UL112 




UL117 


UL114 




UL114 






DNA replication and processing UL44 


UL114 


UL98 


UL89 


UL21A 




UL97 






UL98 






UL114 






Capsid assembly and genome encapsidation UL80 




UL85 


UL85 




UL98 


UL86 






UL97 






UL98 






UL104 






Maturation and release UL4 


UL25 (48 h p 


i.) UL25 (72 h p.i.) 


UL25 


UL26 


UL4 


UL82 


UL43 


UL32 


UL83 


UL50 


UL82 


UL94 


UL55 


UL83 


UL97 


UL69 


UL99 


UL99 


UL96 




UL115 


UL115 






UL128 




Immune evasion RL11 


RL11 




UL119 


UL22A 




UL141 


UL33 




UL144 


UL119 




U373MGs was much less than that in RPE-ls, and the 
burst of intracellular infectivity in both cell lines at 72 h 
p.i. remained largely cell-associated thereafter. Thus, it 
appears that the bottlenecks affecting CRV yield remained 
largely in force for U373MGs, perhaps because the cognate 
functions involved in virion maturation and release were 
provided by viral structural proteins that are needed in 
high abundance. 

It is notable that HCMV genes specifying immune evasion 
functions were differentially regulated in RPE-ls and 
U373MGs (Table 3). Although immune evasion genes 
affect viral pathogenesis in vivo, the effects of down- 
regulation on viral replication in cell culture are unknown. 

Cell type-specific factors that may modulate 
HCMV gene expression 

Cell type -dependent differences in the interactions between 
host and virus are probably responsible for modulation of 
HCMV gene expression. The eye, and to a lesser extent the 
brain, are immune privileged organs, and their response to 
pathogens differs from that of other tissues and organs. 


The RPE layer is part of the blood-retina barrier and 
plays an important part in immune privilege (Holtkamp 
et al, 2001). The immune response of the eye to HCMV 
infection is weak, usually lacks neutrophil involvement, 
and results in persistent infection of retinal tissue (Cinatl 
et al, 2000) and low levels of viral replication (Scholz 
et al, 2003). The observed minimal CPE in RPE-ls and 
maintenance of cell-to-cell contact in the monolayer are 
consistent with preservation of the retinal barrier. Immune 
privilege in the eye may also explain the downregulation oi 
some HCMV immune evasion genes in RPE-ls. 

Binding of the HCMV envelope glycoproteins gB and gH 
to receptors on fibroblasts, and most other cell types, 
causes rapid induction of NF-kB activity, resulting in 
activation of the viral major immediate- early promoter 
(MIEP) and cellular genes that function in cell-mediated 
inflammatory responses, including those involving 
cyclooxygenase-2 (COX-2) and interferons (Eickhoff & 
Cotten, 2005). However, the NF-kB pathway and COX-2 
are not induced in HCMV-infected RPE cells (Cinatl et al, 
2001; Scholz et al, 2003). COX-2 mediates the synthesis oi 
prostaglandins, which are involved in acute and chronic 
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inflammatory and innate immune responses. Thus, avoid- 
ance of NF-kB activation in RPE cells would circumvent 
prostaglandin-associated inflammatory responses in the 
eye and limit viral replication in the retina (Scholz et al, 
2003). Inhibition of COX-2 activity in HCMV- infected 
fibroblasts results in reductions in viral mRNA levels, viral 
DNA synthesis, viral growth and cell-to-cell spread of 
infectivity (Schroer & Shenk, 2008; Zhu et al, 2002). Thus, 
low levels of COX-2 in RPE- Is may contribute to impaired 
viral replication. 

An investigation of HCMV infection in the brain showed 
that neoplastic astroglial cells are permissive for HCMV, 
while non-neoplastic cells are not (Poland et al, 1990). The 
ability of HCMV to grow in neoplastic astrocytomas is 
considered to be dependent upon the degree of cellular 
differentiation, suggesting that differentiation-dependent 
cellular genes may control the HCMV transcriptional 
cascade. HCMV strain Merlin developed a CPE in 
U373MGs that was characterized by large numbers of 
swollen and syncytial cells. Nevertheless, cell-to-cell contact 
appeared to be maintained, suggesting the preservation of a 
barrier cell layer, as in the case of RPE- Is. Infectious viral 
yields from U373MGs were very low compared with those 
from HFFF-2s and RPE- Is, and viral genomic yield was 
129-fold lower than in infected HFFF-2s at 96 h p.i. Lower 
expression levels of genes with roles in viral maturation and 
release of infectious progeny identify these stages as pro- 
bable bottlenecks (Table 3). However, the apparent dereg- 
ulation of the HCMV transcription cascade in U373MGs 
probably accounts for the greater part of the effects on viral 
yields and growth kinetics. 

HCMV utilizes many cellular transcription factors to 
modulate viral gene expression, and cell-type differences 
in the repertoire and levels of these factors may impact 
viral transcriptome activity. The MIEP has multiple 
binding sites for a wide range of factors (Stinski & Meier, 
2007), binding sites for certain factors have been identified 
in UL4 (Chen & Stinski, 2000), UL98 (Adam et al, 1995) 
and UL112 (Rodems et al, 1998), and many other HMCV 
genes are associated with p53-binding sites (Hannemann 
et al, 2009). In HCMV-infected cells, p53 is sequestered in 
viral replication compartments in the nucleus (Fortunato 
& Spector, 1998; Rosenke et al, 2006) and modulates the 
expression of at least one HCMV gene (UL94) (Casavant 
et al, 2006; Wing et al, 1998). The HCMV genome con- 
tains 21 consensus-binding sites for p53, located either in 
the promoters or coding sequences of genes that have roles 
in gene regulation, DNA replication and packaging, and 
maturation and envelopment of viral particles (Rosenke 
et al, 2006). U373MGs have a mutation (R273H) in p53 
(Van Meir et al, 1994). Studies comparing HCMV gene 
expression in fibroblasts either expressing or not expressing 
p53 led to the identification of 22 differentially expressed 
genes, 11 of which were downregulated in non- expressing 
cells, resulting in delayed and reduced infectious viral 
yields (Hannemann et al, 2009). Of the HCMV genes that 
were differentially expressed in RPE-ls or U373MGs, 17 



contain potential p53-binding sites (Hannemann et al, 
2009; Rosenke et al, 2006) and may respond to cell type- 
dependent differences in p53 activity. The unusually rapid 
expression of viral E and L genes in U373MGs may be a 
consequence of the mutated p53 protein. In transient 
transfection studies, the mutated p53 protein activated the 
HCMV MIEP, and the additional presence of CREB 
enhanced this activity. Since p53 mutant-mediated trans- 
activation requires only a functional TATA box (Deb et al, 
1992), it is possible that viral promoters are activated non- 
specifically in U373MGs, regardless of their normal 
temporal kinetics. 

HCMV has evolved the capacity to replicate in various cell 
types that differ in expressed cellular genes and perform 
specialized functions in the context of particular organs 
or tissues. It remains to be determined whether the 
transcriptome activity of strain Merlin in HFFF-2s, RPE- 
ls and U373MGs is mirrored in the cognate cell types in 
vivo. We have shown that expression of particular HCMV 
genes in cultured cells is modulated in response to 
different intracellular environments. Our data suggest that 
the imperatives of immune privilege in RPE-ls and 
U373MGs and differences in expressed cellular factors 
result in reduced transcription from multiple viral genes. 
However, it is highly probable that additional cell type- 
dependent factors not identified here also contribute to 
cell type-specific differences in viral growth kinetics, 
including differential translation of transcripts, post- 
translational processing and intracellular translocation of 
viral products. 

METHODS 

Virus and cell lines. HCMV strain Merlin (Dolan et al, 2004; 
RL13~UL128~) was used. Monolayers of HFFF-2s (ECACC; 
86031405) and U373MGs (ECACC; 89081403) were grown in 
Dulbecco's modified Eagle's medium (DMEM) supplemented with 
10% FCS. hTERT-immortalized RPE-ls (C4000-1; Clontech) were 
grown in DMEM/Ham's F-12 supplement containing 10% FCS. 

Susceptibility to infection. HFFF-2s, RPE-ls or U373MGs (2 x 10 5 
cells) were infected with HCMV at 6 p.f.u. per cell, and the number of 
infected cells in two fields of view (nominally 500-600 cells per field) 
was estimated by immunofluorescence at 48 h p.i. by using an UL44 
antibody (CCH2-FITC; Dako). 

Determination of viral genome yields. HFFF-2s, RPE-ls and 
U373MGs (4 x 10 5 cells) were infected with HCMV at 6 p.f.u. per 
cell, and whole cell DNA was harvested at 24, 48, 72 and 96 h p.i. 
DNA was also harvested at 48 h from MI cells. The monolayers were 
trypsinized and cells were pelleted at 550 g for 10 min at room 
temperature. Infected cell DNA extracts were prepared by using a 
FlexiGene kit (Qiagen). Viral genome yield was estimated by qPCR 
using UL130 gene primers ( 5 ' - GCGAGGGATAGAGAAAAGG- 
ACAG-3' and 5 ' -CCGTGGTCGACGCTAACAG-3 ' ) and a TaqMan 
probe (5 ' -6-FAM-CGGTTTGGAATACGTCAGT-MGB-3 ' ) . Each 
experiment involved amplification of MI, infected cell and control 
plasmid DNAs, with four reactions per sample. The plasmid, which 
contained the UL130 ORF, was used to generate a standard curve for 
determination of HCMV genome copy number in infected cell 
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samples. Samples were amplified (40 cycles at 95 °C for 3 s followed 
by 60 °C for 30 s) in an Applied Biosystems 7500 Fast Real-time 
PCR System, and data were analysed using proprietary software. 

Preparation of RNA from Ml and HCMV-infected cells. For 

microarray analysis, HFFF-2s, RPE-ls and U373MGs were infected 
simultaneously with a single HCMV preparation, and RNA was 
harvested at 12, 24, 48 and 72 h p.i. For each cell line, 12 monolayers 
(each 6 x 10 6 cells) were infected at 6 p.f.u. per cell, thus providing 
three biological replicates for each time point. MI cultures were 
harvested at 72 h p.i. Cells were solubilized in lysis buffer containing 
1 % 2-mercaptoethanol, and total cell RNA was extracted using an 
RNeasy kit (Qiagen). After treatment with DNase I (Invitrogen), the 
quality and purity of RNA preparations were confirmed by lack of 
smearing of the 28S and 18S rRNA bands after electrophoresis on 1 % 
agarose/2.2 M formaldehyde gels. Elimination of DNA from the 
preparations was confirmed by failure to amplify a cellular DNA 
sequence (from the lactate dehydrogenase gene) by PCR. 

Northern blot analysis. RNA (10 ug per track) was loaded onto a 
1 % agarose/2.2 M formaldehyde gel, separated by electrophoresis, 
transferred by blotting to Hybond-N+ membrane (Amersham 
Biosciences), and fixed by irradiation in a UV-cross-linker 
(Stratagene) at 12 000 J cm~ 2 . DNA probes were prepared from 
PCR-generated templates by using a Rediprime II random prime 
labelling kit (Amersham Biosciences), with direct incorporation of 
[a- 32 P]dCTP (50 uCi) (Amersham Biosciences), diluted and de- 
natured in Rapid-hyb buffer (Amersham Biosciences), and hybridized 
to the membrane overnight at 68 °C. The hybridized membranes 
were washed and labelled RNA bands were visualized by using a 
Bio-Rad Personal FX phosphorimager. 

Microarray probe design and fabrication. Specific 60-mer 
oligonucleotide probes (Table S2) were generated for HCMV ORFs 
(GenBank accession no. AY446894.2) according to the protocol 
published previously (Ebrahimi et al., 2003), synthesized by MWG- 
Biotech (USA), and printed in triplicate on Corning GAPS II 
microscope slides, with 16 copies of the microarray on each slide. 
Probes for two Bacillus subtilis genes (spoOB, GenBank accession no. 
M24537, and lys, GenBank accession no. X17013) were utilized in both 
sense and antisense orientations in order to assess array specificity and 
sensitivity and to serve as normalization controls. The array was 
validated by hybridization against MI and infected cell cDNAs at a 
range of temperatures in the presence or absence of spiked controls. 

Synthesis of cDNA for microarray hybridization. RNA aliquots 
(25 ug) were spiked with B. subtilis control mRNAs (50 ng each), mixed 
with 1 ug random hexamer primers (Invitrogen), incubated at 70 "C 
for 10 min, and cooled on ice. Labelled cDNAs were then produced by 
incubating the samples at 42 °C for 2 h with Cy3-dCTP (Roche), 
dNTPs (Roche) and Superscript II reverse transcriptase (Invitrogen). 

Microarray hybridization, data collection and data processing. 

Labelled cDNAs were prepared for hybridization as described 
previously (Ebrahimi et al, 2003). Hybridization was performed in 
a humidified chamber (Genetix). Array specificity was confirmed by 
hybridization with cDNAs prepared from MI- or HCMV-infected cell 
RNAs. Arrays were scanned at 543 nm with a pixel resolution of 
10 um using ScanArray Express hardware and associated software 
(Perkin-Elmer). The background intensity value was subtracted from 
the corresponding spot intensity, and the mean value for each probe 
was calculated («=3). The linear dynamic range was determined by 
using published methods (Ebrahimi et al., 2003). 

The signal threshold that gave a false positive rate of 5% was 
determined from receiver operating characteristic analysis of the 
control data. Datasets were then scaled against the 75th percentile of the 



positive-control signal. Probes with signal values falling below the 
signal threshold cut-off point were assigned a signal value of 0.01 to 
facilitate log-transformation and further analysis of data. ORFs UL48A, 
US1 and US34A gave signal values below the detection limits of the 
assay for all time points in each cell line, and were not analysed further. 

Array variability was checked by measuring the median signal 
intensity for a given cell line and time point p.i. on log-transformed 
data. Only arrays where total median signal intensity values were 
within threefold of each other were selected for further analysis. To 
identify differentially expressed genes, mean log 2 signal intensity 
values were calculated for all viral ORFs at each time point in the 
three cell lines. Where values were >0.01 in both the control (HFFF- 
2s) and test cells (RPE-ls and U373MGs), signal intensity fold change 
(FC) values were computed. The probability of cell type-dependent 
differential gene expression (P<0.05) was then determined by 
Student's t-tests of FC values for each time point (Table SI). 
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